Universality of Many-Body States in Rotating Bose and Fermi Systems 
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We propose a universal transformation from a many-boson state to a corresponding many-fermion 
state in the lowest Landau level approximation of rotating many-body systems, inspired by the 
Laughlin wave function and by the Jain composite-fermion construction. We employ the exact- 
diagonalization technique for finding the many-body states. The overlap between the transformed 
boson ground state and the true fermion ground state is calculated in order to measure the quality 
of the transformation. For very small and high angular momenta, the overlap is typically above 
90%. For intermediate angular momenta, mixing between states complicates the picture and leads 
to small ground-state overlaps at some angular momenta. 
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The properties of rotating many-body systems have 
been the topic of intense study, theoretical as well as 
experimental ^3M Jii 11.0, i S [M EH, M El, 
E3, EM EB, E3,lll, E3, Ho]- particular, Bose-Einstein 
condensates have been of great interest since the advent 
of the laser-cooling technique [2l|. It is by now a well- 
established fact that quantum many-body systems un- 
der rotation form quantized vortices, a property that has 
been known in the contexts of superconductivity and su- 
perfluidity [22, HI]. In earlier works, the formation of vor- 
tices in few-body systems was studied quite extensively, 
including both boson and fermion systems , as 

well as their two-component generalizations ^24 . i25i . i26] . 
In Ref . (I we noted that some properties of vortex forma- 
tion in few-body systems of repulsively interacting par- 
ticles in the lowest Landau level are universal, not only 
with respect to the details of the interaction (the boson 
system was studied both with the long-range Coulomb 
interaction and with a short-range, contact interaction), 
but also with respect to the statistics of the constituent 
particles. We found remarkable similarities between the 
yrast spectra of the boson and fermion systems, and we 
also noted that vortices enter at very specific, and cor- 
responding, angular momenta in the boson and fermion 
systems, respectively. 

In the present paper, we formalize this boson-fermion 
universality by means of a transformation from a bosonic 
many-body wave function to a corresponding fermionic 
wave function. The transformation is inspired by the 
Laughlin wave function [13], and forms a direct paral- 
lel to the Chern-Simons transformation in the theory 
of the fractional quantum Hall effect [2^, |2^ and Jain's 
composite-fermion picture |3(^] . (The latter was recently 
successfully applied to the problem of small quantum 
dots in magnetic fields, and other small quantum systems 
at high angular momentum [1, E3, El El El, El, ES, iS] ■ ) 
This allows us to investigate quantitatively to what ex- 
tent this universality holds in a general comparison be- 
tween few-body boson and fermion systems in the lowest 
Landau level. The particles may be electrons in a quan- 
tum dot (in the case of fermions) or optically or mag- 



netically trapped ions at ultra- low temperatures. In this 
work we concentrate on harmonically confined. Coulomb- 
interacting particles in two dimensions. 

This paper is organized as follows: The next section 
gives an introduction to the mathematics of the pro- 
posed boson-fermion universality of many-body systems 
in the lowest Landau level. This section also introduces 
the concepts used throughout the paper, and defines 
a mathematical transformation of a many-boson wave 
function into a corresponding many-fermion wave func- 
tion. Section |TT] gives a brief description of the exact- 
diagonalization method in the lowest Landau level, and 
of the implementation of the boson-fermion transforma- 
tion. The reader who is familiar with exact diagonaliza- 
tion may skip directly to section IIIIl where our results 
are presented and discussed. Some concluding remarks 
are given in section [TVl 



I. INTRODUCTION TO BOSON-FERMION 
UNIVERSALITY 

In order to formulate a mathematical expression for 
the direct comparison between boson and fermion many- 
body wave functions, we use complex coordinates for the 
two-dimensional plane, z — x + iy. The general wave 
function in the lowest Landau level of one particle in a 
harmonic-oscillator potential is 4'e{z) oc z^e~'^' , where i 
is the angular momentum. (Atomic units and one fixed 
frequency lu — 1 of the harmonic confinement are used 
throughout the paper.) In the boson system at zero an- 
gular momentum, all particles reside in the £ = state, 
and the many-body wave function is 



(1) 



As the system is set rotating, particles are lifted from 
the £ — single-particle state in order to carry angular 
momentum. A single vortex is formed at the center of 
mass when the total angular momentum of the system 
equals the number of particles, in which case all particles 
are in the £ = 1 state. As angular momentum increases 
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further, more vortices successively enter the system. As 
detailed in Ref. 3, a vortex-generating state carrying 
n vortices can be obtained from Eq. ([T]) by successive 
multiplications with symmetric polynomials: 
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This wave function describes a state with n vortices 
evenly spaced on a ring with radius a centered at the 
origin, and may be used to obtain a trial many-body 
wave function [8| . 

The same line of reasoning can be repeated completely 
analogously for fermions. However, in the lowest Landau 
level the angular momentum of a many-fermion system 
cannot be zero. Instead the smallest possible total angu- 
lar momentum is achieved by putting one fermion in each 
of the N lowest single-particle states with single-particle 
angular momenta ranging from to — 1. This state 
is called the maximum density droplet (MDD) and is the 
fermion equivalent of the zero-angular-momentum states 
for bosons. The wave function of the MDD is given by 
the Laughlin wave function with filling factor one [sij : 
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In the construction of the vortex-generating state for the 
fermion system, the boson condensate is replaced by the 
MDD. 

Now we compare Eq. ([Ij with Eq. ([3]) and note that 
the latter can be obtained from the former simply by 
multiplication with the polynomial 
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Making the assumption that the transformation between 
the lowest angular momentum states, V^mod = ^^^a i 
holds whenever is replaced by ^mddi arrive at 
a very general transformation from any bosonic many- 
body wave function to a corresponding fermionic wave 
function, shifted in angular momentum by exactly Lmdd: 



(5) 



The degree of universality of the structure of the many- 
body wave function between boson and fermion systems 
will be reflected in how well the transformation ^ re- 
produces the true fermion wave function. The straight- 
forward way of determining this is to calculate the wave 
function of a given many-fermion system both directly, in 
order to obtain the true wave function, and by transform- 
ing the corresponding boson state according to Eq. 



The overlap between the true wave function and the 
transformed boson wave function can the be calculated: 



\\ iJMDD+i^B I J^b/\ 
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The argument given above for the boson-fermion cor- 
respondence is entirely heuristic: it is based on observing 
the common features of fermion and boson wave func- 
tions in the lowest Landau level. The correspondence 
formalized in Eq. ([5]) is ultimately justified by the ex- 
plicit demonstration of its performance in section IIIII 
However, it is reasonable to anticipate the existence of a 
boson-fermion transformation of this type from boson- 
Chern-Simons theory [2^ [2^ and composite-fermion 
theory (29l . [30| . A fermionization of repulsively interact- 
ing bosons in the lowest Landau level has been described 
by several authors 0, UIl, [H, fisl . [TgI . Iit\ . In particular, 
Cazalilla et al. , Chang et al. as well as Regnault 
et al. [iBl explicitly map interacting bosons onto non- 
interacting, spinless fermions by means of a composite- 
fermion construction. It is known that the mean-field ap- 
proximation of non-interacting composite fermions loses 
accuracy as angular momentum increases, and effects of 
residual interactions become important [§, [ll|, [H, . In 
the case of a harmonic interaction, the transformation ([5]) 
between interacting bosons and interacting fermions was 
rigorously derived recently by Ruuska and Manninen [3^ . 

In the Chern-Simons approach to the fractional quan- 
tum Hall effect, the many-body wave function of the 
electrons of the two-dimensional electron liquid is trans- 
formed into a wave function of particles moving in 
an effective magnetic field. The transformation is defined 
as follows 1291: 
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The prefactor determines the number of flux quanta as- 
sociated with each particle, and it is symmetric or anti- 
symmetric with respect to particle interchange depend- 
ing on whether p is chosen even or odd. In the sim- 
plest possible model, we choose p = 1. Since the original 
electron wave function is fermionic, this means that the 
Chern-Simons wave function describes a corresponding 
boson system. The transformation between the boson 
system and the fermion system is the given by the factor 

Jain took the Chern-Simons approach further to a 
more sophisticated ansatz in the lowest Landau level by 
introducing his composite-fermion theory (29l . [30| . In 
this approach, the power p of the Chern-Simons factor 
is explicitly taken to be an even integer, thereby mak- 
ing the composite particles fermions. Also, the prefactor 
is simply taken to be the Jastrow factor W^^jizi — Zj)^, 
which attaches vortices rather than flux tubes to the par- 
ticles [1^. The wave function of the original fermions is 
retained by projection to the lowest Landau level: 



Vl[{z,~z,f^CF, 



(8) 
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where V is the projection operator into the lowest Landau 
leveL 

Jain's construction and the bosonic Chern-Simons 
wave function immediately suggest a mapping between a 
boson system and an equivalent fermion system at some 
other magnetic field (or, equivalently, angular momen- 
tum) of the form described by Eq. For a harmonic 
particle-particle interaction, this relation is exact, and 
can be derived analytically by constructing a spectrum- 
generating algebra, and applying the corresponding lad- 
der operators to the ground-state wave functions ^ b = 1 
and — Y\.i<:j {^i ~ ) respectively [32| . 

The numerical study presented in this paper shows 
that the transformation is very general and yields a good 
approximation for the fermion wave function in the low- 
est Landau level, also for Coulomb-interacting particles. 



Boson wave function Transform 

0.685 1 04000) -5= 0.895 1 01111000) 

-0.168 1 30001) -0.400 1 10110100) 

H-0.194|21010) /^^^^^^^ +0.138 1 11001100) 

+^.m\\2<d2m)^^^-^'^'''^^^ +0.118 1 11010010) 

-0.559 1 12100) -0.0775 1 11100001) 

FIG. 1; As an example, we show the transformation of the 
Lb ~ 4 boson wave function into a fermionic wave function 
with Lf = 10. One term in the boson wave function may 
contribute to the weight of several terms in the transformed 
wave function. Correspondingly, terms in the transformed 
wave function may receive contributions from more than one 
term in the original wave function. 



II. THE EXACT-DIAGONALIZATION 
METHOD AND THE BOSON-FERMION 
TRANSFORMATION 

In order to calculate the overlap, Eq. the many- 
body wave functions must be obtained with sufficient ac- 
curacy. In this study we take as our model system har- 
monically confined. Coulomb-interacting particles, and 
we solve both the fermion and the boson problems with 
the exact-diagonalization method. The boson wave func- 
tion is transformed into a corresponding fermionic many- 
body wave function by means of Eq. ([5]), which we have 
implemented numerically to operate on a bosonic wave 
function in the Fock-Darwin basis. 

Since the model system used in this study is that 
of a harmonic confining potential, a natural choice of 
single-particle basis is the set of eigenfunctions of the 
two-dimensional harmonic oscillator in the lowest Lan- 
dau level: 
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The many-body wave function is then built in the Fock- 
Darwin basis formed from these single-particle states, 
and the Hamiltonian (up to an additive constant) is 
straightforwardly written as 



i i.j,k,l 

where the interaction matrix element is 



(10) 



Ut,j,k,i = (ij ^(r,r') 



kl 



{vWj {r')U{r, r')yjk (r)V;, (r') dr dr' . 

(11) 

Since we are working with many-body systems at given 
total angular momentum L, the one-body term yields 
the same energy contribution for all allowed Fock states. 
Therefore we may drop this term, and diagonalize only 



the interaction part of the Hamiltonian. Working in the 
lowest Landau level already guarantees that we have a 
finite number of possible Fock states for any given to- 
tal angular momentum, and in this study we have used 
all of these in the basis. Thus no further truncation of 
the Hilbert space is done. The resulting matrix is di- 
agonalized numerically, yielding the interaction-energy 
spectrum and the many-body eigenstates. The diago- 
nalization is performed using the Lanczos algorithm [33| , 
except at the smallest angular momenta, where a full di- 
agonalization is needed due to the small matrix sizes. 

The exact-diagonalization method outlined above is 
applicable to both fermions and bosons; only the commu- 
tation relations of the creation and annihilation operators 
differ, and the many-body Fock-Darwin basis states must 
be given the appropriate symmetry properties (Slater de- 
terminants for fermions, product states for bosons). 

The numerical implementation of the transformation, 
Eq. ([5]), requires some care. A general bosonic many- 
body wave function consists of a linear combination of a 
(possibly very large) number of Fock states. Upon multi- 
plication by the determinant (Eq. each bosonic 
Fock state is transformed into a fermionic many-body 
wave function, which is in general not a basis state in the 
fermion Fock-Darwin basis, but instead must be repre- 
sented as a linear combination of fermionic basis states. 
Thus, each Fock state in the boson wave function will 
yield several Fock states in the final fermionic wave func- 
tion, and correspondingly, the weight of each Fock state 
in the final fermionic wave function may receive contri- 
butions from several Fock states in the original bosonic 
wave function. 

For example, the four-particle bosonic state 1 202) 
(with angular momentum Lb — ^) transforms into the 
fermionic wave function 0.79 1 110011) - 0.56 1 101101) + 
0.25 I 011110), a linear combination of three states from 
the Fock basis (total angular momentum Lp = 10). The 
true wave function for the four-boson system with angu- 
lar momentum Lb = 4 is a linear combination of five 
basis states, as depicted in Fig. [TJ The state in our 
toy example appears in the fourth term. Under appli- 
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cation of the transformation, Eq. ([5]), this wave func- 
tion transforms into the fermionic wave function shown 
in the right column of Fig. [1] As indicated by arrows, the 
bosonic state | 202) contributes to the first three terms of 
the transformed wave function. Contributions come also 
from other terms in the original wave function. For ex- 
ample, contributions to the first term come from all terms 
of the bosonic wave function. The contribution from the 
boson basis state | 04000) is also indicated in the figure 
(in fact, this particular term in the boson wave func- 
tion contributes only to the first term of the transformed 
wave function). This example shows that it is necessary 
to carefully collect the weights of the basis states of the 
fermionic wave function resulting from the transforma- 
tion. 



III. RESULTS 

We perform calculations of energies, wave functions 
and ground-state overlaps systematically over a wide 
range of angular momenta for systems containing be- 
tween four and eight Coulomb-interacting particles. 

Fig. [5] shows the overlap of the ground state of a four- 
fermion system with the transformed ground state of a 
corresponding four-boson system, for total fermion an- 
gular momenta in the range Lp = 18-62. In the small 
four-particle system, the overlap is very good, close to 
100% for all angular momenta. The top and middle pan- 
els of the figure show the yrast spectra of the five lowest 
states for the boson and fermion systems respectively. 
The similarity of the yrast lines is readily apparent from 
this figure. Both systems show a four-fold periodicity, 
and the kinks in the yrast line — where there is also a 
relatively large gap between the ground state and the 
first excited state — appear at corresponding angular mo- 
menta, Lp = Lb + Lmdd- The overlap oscillates with 
the same four-fold period for high angular momenta, with 
the overlap having a local maximum whenever there is a 
kink in the yrast line. 

Looking at a system with more particles, the picture 
becomes more complicated. Fig. |3] shows the ground- 
state overlap in the five- and six-particle systems over 
the entire range of total angular momenta studied, Lp — 
Lmdd = 10 to L^? = 85 for five particles and Lp ~ 
Lmdd — 15 to Lp — 89 for six. At the smallest angular 
momenta, the overlap is always 100% due to the fact that 
the Hilbert space is so small that only one or a few wave 
functions are possible. It is interesting to note that also 
in the limit of high angular momentum, the overlap be- 
tween the true fermion ground state and the transformed 
boson ground state tends towards 100%. This shows that 
the performance of Eq. ([5]) is not a consequence of the 
Laughlin wave function being a good approximation. In 
fact, the overlap between the Laughlin wave function and 
the exact wave function decreases with decreasing filling 
factor 

Fig. [3] also shows some cases where the overlap drops 



to zero. These cases will be discussed in some detail later 
in this section 

Away from the extremes of very small and very large 
angular momenta, the overlap drops and begins to fluc- 
tuate. In the case of five particles (top panel of Fig. [3]) 
this is seen as a general decrease of the overlap for angu- 
lar momenta in the range Lp = 14-35 together with two 
drops to zero a,t Lp = 16 and Lp = 46. This pattern is 
repeated in a much more complicated picture in the six- 
particle case, shown in Fig.[31 bottom panel. Again, there 
is a general decrease in the overlap as the total angular 
momentum increases beyond the small angular momenta 
where the overlap is large due to the small Hilbert space. 
As the angular momentum increases beyond Lj? « 40 
there is a region where the overlap tends to be either 
large or vanishing. At the very largest angular momenta 
studied, the overlap settles at over 90%. 

The behavior described for the five- and six-particle 
systems is echoed in smaller and larger systems as well. 
In the four-particle system (Fig. (2) the overlap is well 
above 90% for all angular momenta studied, Lp = 6-82. 
However, for angular momenta in the range Lp = 10- 
40 the overlap oscillates between 96% and 99%, and for 
Lp = 24 dropping to 93%, before returning to values 
very close to 100% for larger angular momenta. In sys- 
tems with seven and eight particles (Fig. d]) the overlap 
oscillates wildly over most of the studied region of an- 
gular momentum, where overlaps at or above 80% are 
interspersed with drops to about 50% or below. There is 
a tendency towards better overlaps in the high angular 
momentum limit of the data. 

We may conclude from these observations that the 
transformation defined in Eq. ([S]) from a bosonic many- 
body wave function to a fermionic wave function for 
the same number of particles and angular momentum 
Lp = Lb + ^MDD performs well in the low and high an- 
gular momentum limits, where good performance may be 
expected due to the restricted Hilbert space and particle 
localization respectively. For intermediate values of the 
angular momentum, the performance of the transforma- 
tion needs a deeper study. 

The poor performance of Eq. ([5]) for certain combi- 
nations of particle number and angular momentum is a 
result of a restructuring of levels in the yrast spectrum 
between the boson and fermion systems. The result is 
that the boson ground state will have a different structure 
than the fermion state, and therefore so will the trans- 
formed boson wave function. In order to study the par- 
ticle configurations corresponding to the different many- 
body wave functions, wc look at the pair- correlation func- 
tion 



7/'t(r')v;t(r)v3(r')7/'(r 



(12) 



which gives the probability of finding a particle at r, given 
that there is a particle at r'. 

We may compare the pair-correlation functions ob- 
tained from the true fermion wave function and from the 
transformed boson wave function. Six particles at total 
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FIG. 2: Yrast spectra and overlaps for a four-particle system. (A) Yrast spectrum of four harmonically confined bosons. The 
five lowest states for each angular momentum are included. (B) The corresponding yrast spectrum for four fermions. (C) 
Overlap of the transformed boson ground state with the fermion ground state, calculated from Eq. ((6)1. 



fermionic angular momentum Lp = 57 and Li? = 58 are 
two examples from the intermediate angular momentum 
range where the transformed boson ground state some- 
times reproduces the fermion ground state very well, but 
sometimes fails spectacularly. The former is the case at 
Lp = 57, while the latter is true for Lp = 58. Fig. [5] 
shows the pair-correlation functions of the boson wave 
function (top panel), the transformed boson wave func- 
tion (middle panel) , and the fermion wave function (bot- 
tom panel) for Lp — 57 (left column) and Lp — 58 (right 
column). The reference particle in each plot is indicated 
by a black dot in the figure. Its position is chosen such 
that it sits at the distance rmax from the origin where the 
particle density has its maximum. 

At Lp = 57 the particle structure is extremely well 
reproduced by the transformation from the bosonic state 
with Lb = Lp—LuBB = 42, and the two pair-correlation 
functions from the true fermion wave function and the 
transformed boson wave function are almost indistin- 
guishable. The true fermionic many-body wave function 
displays an internal structure of localized particles, where 



the six particles of the system are equidistantly situated 
on a circle. This is one of two classically (meta-)stable 
configurations for six coulomb-interacting particles in 
a harmonic confining potential \S5i\ . This structure is 
also seen in the the bosonic wave function, and the ef- 
fect of multiplication with , apart from changing the 
particle-exchange symmetry, is to push the particles out- 
ward. 

The other classically stable configuration of six equal 
electrical charges is a ring of five particles with the sixth 
particle sitting precisely at the center of the ring. In 
the classical system, this configuration actually has a 
slightly lower energy than the ring of six particles [1^. 
These two configurations, (6,0) and (5, 1), are known to 
compete in quantum-mechanical six-body systems with 
repulsive coulomb interaction (sgI . [37j . In an electron 
system where the spin degree of freedom is not frozen 
out, the (6,0) configuration tends to be favored (unless 
the system is very dilute), because of frustration in the 
(5, 1) configuration [H, [13, Ull . In the present study, we 
deal with systems of spinless bosons and of spin-polarized 
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FIG. 3: Overlap between the true fermion ground state and 
the transformed boson ground state as a function of total 
angular momentum for five (top panel) and six (bottom panel) 
particles. 
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FIG. 4: Overlap between the true fermion ground state and 
the transformed boson ground state as a function of total 
angular momentum for seven (top panel) and eight (bottom 
panel) particles. 



fermions. Therefore, spin frustration does not come into 
play, and the two classical configurations compete. When 
the angular momentum is increased just by one unit to 
Lp = 58, the fermion ground state is the (5, 1) configura- 
tion, as shown in the bottom right panel of Fig. \5\ In the 
boson system, however, another change takes place with 
the corresponding increase of one unit of angular mo- 
mentum: the particles remain seated on one ring of six 
particles, but the degree of localization decreases. Corre- 
spondingly, the multiplication with expands the state 
to a ring of six slightly smeared out maxima, a state that 
is orthogonal to the true ground state. 

Indeed the orthogonality of the transformed state to 
the true fermion ground state is no coincidence, as in fact 
the transformed bosonic ground state corresponds to the 
first excited state of the fermion system, and correspond- 
ingly the fermion ground state is well reproduced by ap- 
plying the transformation to the first excited state of the 
boson system. Table U shows the overlap of the fermion 
ground-state wave function with each of the transformed 
wave functions of the five lowest boson states. The same 
table also shows the overlap of the transformed boson 
ground-state wave function with each of the five lowest 
states of the fermion system. As can be seen from the 
table, the transformed bosonic ground state almost per- 
fectly reproduces the first excited state of the fermion 
system. Plotting the pair-correlation functions, as done 
in Fig. [SI confirms this picture. This figure shows the 
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FIG. 5; Pair-correlation functions for the six-particle sys- 
tem at total angular momentum Lp ~ 57, corresponding to 
Lb ~ 42 in the boson system (top row), and ai Lp —58 
(bottom row). The panels in each row show, from left to 
right, the pair-correlation function of the boson ground state, 
the transformed boson ground state, and the fermion ground 
state, respectively. The black dot indicates the position of the 
fixed reference particle 

pair-correlation functions of the first excited state in the 
boson system and the transformed version of this state, 
and the first excited fermionic state. These plots may 
be compared with the corresponding ground-state pair 
correlations in Fig. [51 

N = 6, Lp = 58 discussed above is a particularly clean 
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TABLE I; Overlap of the fermion ground state with the trans- 
formed wave function of each of the five lowest states of the 
boson system (middle column), in the system with N — G 
particles and fermionic angular momentum Lf — 58. The 
table also shows the overlap between the transformed boson 
ground state with each of the five lowest fermion states in the 
same system (right column) 
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FIG. 6: Pair-correlation functions for the first excited state in 
the six-particle system at total angular momentum Lp ~ 58, 
corresponding to Lb = 43 in the boson system. The panels 
show the pair-correlation function of the boson wave function 
(left), the transformed boson wave function (center), and the 
fermion wave function (right), respectively. The black dot 
indicates the position of the fixed reference particle. 



example where the many-body configurations exchange 
places between the fermion and the boson spectra. More 
generally, states may mix, such that one fermionic state 
is reproduced by some linear combination of transformed 
boson states. Such mixing of states between the bosonic 
and fermionic systems accounts for the cases where the 
overlap between the ground states is neither close to 
100% nor vanishing. One example is TV = 6 particles at 
angular momentum Lp = 53, where the overlap between 
the ground states is 68.32%. In this case the fermionic 
ground state has a considerable overlap with all of the 
three lowest bosonic states. This is shown in Table |lll 
The pair-correlation functions are shown in Fig. [71 The 
transformed boson ground state displays a smeared-out 
(5, 1) configuration, while the fermion ground state is a 
(6, 0) configuration distorted such that there is also an 
increased particle density at the origin. 

The described mixing of the order of the states between 
bosonic and fermionic system is seen also in the other 
cases were the overlap between the ground states is bad, 
also for other particle numbers than six. There are many 
examples at all particle numbers larger than four. One 
particularly drastic example for eight particles is Lp ~ 
57, where the fermion ground state has 91.20% overlap 
with the tranformed wave function of the first excited 



TABLE IL Overlap of the fermion ground state with the 
transformed wave function of each of the five lowest states of 
the boson system (middle column) , in the system with TV = 6 
particles and fermionic angular momentum Lp — 53. The 
table also shows the overlap between the transformed boson 
ground state with each of the five lowest fermion states in the 
same system (right column) 
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68.32 


68.32 
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25.74 


20.71 


2nd 


4.67 


9.29 


3rd 


0.00 


0.00 


4th 


0.00 


0.00 




FIG. 7: Pair-correlation functions for the six-particle system 
at total angular momentum Lf ~ 53, corresponding to Lb = 
38 in the boson system. The panels show the pair-correlation 
function of the boson ground state (left), the transformed 
boson ground state (center), and the fermion ground state 
(right), respectively. The black dot indicates the position of 
the fixed reference particle. 



boson state. However, the transformed bosonic ground 
state does not have any appreciable overlap (1.94%) with 
the first excited fermion state. Instead it corresponds 
more closely to the second excited state (with an overlap 
of 89.91%). 

We have seen that when the particle number increases 
the overlap between the fermion and boson states is still 
good, but in many cases the lowest energy state of the 
fermion system corresponds to an excited state of the 
boson system and vice versa. With increasing particle 
number the exchange of the order of levels becomes more 
important. This is already known in connection with the 
study of vortices in boson and fermion systems [l^ . For 
example, in the boson case the single vortex reaches the 
origin, while in the fermion case the state with one vortex 
at the origin becomes an excited state when the particle 
number is about 10 or greater. This is also true for states 
with multiple vortices: the same vortex states are there 
in both systems, but they appear as the lowest state at 
different (relative) angular momenta. As an example, 
we may study three vortices in a rotating system of 20 
particles. For bosons the angular momentum is 48 and 
the state is the lowest-energy state, while for fermions it 
is the fourth state at angular momentum 238. Figure [8] 
shows the hole-hole correlation |39[ (i.e. pair correlation 
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of vortices) which can be calculated for the fermionic 
wave functions (for the boson system, the transforma- 
tion, Eq. \5\ is first applied to the wave function) . With 
this large number of particles we were not able to per- 
form the computation with the full Fock-Darwin basis. 
The resulting overlap between the converted boson and 
the fermion wave functions is only 88%, but yet the re- 
sulting pair-correlation functions in Fig. [5] are nearly in- 
distinguishable. 




20 FERMIONS, L = 238 



20 BOSONS, L = 48 



FIG. 8: Hole-hole pair-correlation function describing the 
similar vortex localization in boson and fermion systems. The 
boson wave function was first converted to a fermion wave 
function as described in the text. 



pect large overlaps regardless of the type of repulsive in- 
teraction, merely due to the small Hilbert space. It is 
interesting to note, that not only are the overlaps large 
(over 75%), but the overlaps obtained from (5-interacting 
bosons follow those obtained from Coulomb-interacting 
bosons closely. 
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FIG. 9; Overlap between the ground state of six Coulomb- 
interacting fermions and the transformed ground state of six 
bosons with (5-interaction, plotted as a function of total an- 
gular momentum in the fermion system. 



One may speculate that the large overlaps between 
fermionic wave function and transformed bosonic wave 
function occur when the wave function is particularly 
simple. One possible measure of the complexity of the 
fermion wave function within the given single-particle ba- 
sis (here the harmonic-oscillator eigenstates) is the num- 
ber of Fock states needed to make up 50% of its norm. 
For six particles, this number is five or below for angu- 
lar momenta up to Lp = 30, and the overlap is 70% 
and above, but we see very large overlaps also when the 
number of Fock states needed for 50% of the norm of 
the fermion wave function is larger than 40. This shows 
that the good overlaps between fermion state and trans- 
formed boson state are not merely the result of simple 
wave functions, but that the transformation defined by 
Eq. (O can handle also complex many body states. 

The results detailed this far may to some extent be 
robust against changes in the details of the repulsive in- 
teraction between the constituent particles of the system. 
This is suggested by results where we calculate the over- 
lap between the ground-state wave function for Coulomb- 
interacting fermions and the transformed ground-state 
wave function of bosons with a repulsive 5-type interac- 
tion. This is done for the case of six particles at an- 
gular momenta between Lp = 15 and Lp = 45. The 
overlap is plotted as a function of Lp in Fig. [SI This 
figure should be compared with the same range of angu- 
lar momentum in the bottom panel of Fig. [31 For the 
lower half of the range of angular momenta, the overlap 
where (5-interaction is used for the bosons does not de- 
viate substantially from the overlap where both systems 
have Coulomb interaction. For these angular momenta, 
the greater part of the many-body wave function comes 
from only a few Fock states and therefore we may ex- 



These results should be expected from what we know 
about the Laughlin wave function. The Laughlin wave 
function is a good approximation to the exact many- 
body state as long as the angular momentum is not too 
large. Laughlin's construction makes no explicit refer- 
ence to the details of the repulsive interaction between 
the constituent particles [27] . 

For angular momenta Lp > 30 the overlap drops to 
zero except for certain angular momenta, where the over- 
lap is very good (95% and above), and very close to the 
overlap obtained when Coulomb interaction is used. In 
many cases the zero overlap is again simply the cause of 
different orderings of the states in the the two systems, 
but in some cases also a substantial mixing of states was 
found. 



IV. CONCLUSION 

We conclude that for particles in the lowest Landau 
level there is a far-reaching universality between bosons 
and fermions in the properties of the rotating systems. 
This universality may be formulated mathematically as a 
transformation (Eq. ([S])) from a bosonic many-body wave 
function to a fermionic one. These two wave functions 
will differ in total angular momentum exactly by Lmdd, 
the smallest possible angular momentum in the fermion 
system. The transformation produces a very good corre- 
spondence (as measured by calculating the overlap inte- 
gral) between the bosonic and fermionic states when the 
number of particles is small, when the angular momen- 
tum is very small (due to the restricted Hilbert space), 
and when the angular momentum is large (due to local- 
ization into states well described by Laughlin wave func- 
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tions). Away from these extremes, the correspondence 
between boson and fermion states is more compUcated. 
Apart from a general decrease of the the overlaps due 
to the difficulty of reproducing the details of the more 
complicated wave functions as the particle number in- 
creases, there is also fluctuations between large and very 
small overlaps between the ground states. This is due to 
a reshuffling of the levels in the spectrum between the 
boson and fermion systems. In these cases it is not the 
ground state but a low-lying excited state that transforms 



to the fermion ground states. 
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